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Abstract. Double ionization processes triggered by intense linearly polarized laser fields have revealed the paramount role of 
electron-electron collisions. It has been shown [1] that the correlated dynamics of such electronic collisions can be modeled 
by a two dimensional symplectic map which captures the bare essentials of the energy exchange processes occurring during 
a recollision between the two electrons. We investigate linear and nonlinear properties of this map and connect them to the 
<**] ' dynamical organization of phase space and related statistical data. 
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ctf : INTRODUCTION 

Atoms and molecules subjected to intense linearly polarized laser fields have demonstrated the dramatic impact of 
electron-electron correlation in the double ionization process. The ionization channel by which the two electrons 
collide before ionizing (called nonsequential double ionization) has been shown to be the dominant channel for a 
significant range of laser intensity; indeed the corresponding double ionization yields were found to be several orders 
of magnitude higher than the sequential mechanism predicts [3, 4, 5, 6, 7, 8, 9, 10, 11]. The recollision [12, 13] in 
which an electron is first ionized and then hurled back to the core upon sign reversal of the laser field and collides with 
the other electron, has been identified as the "keystone of strong field physics" [2]. It has been confronted successfully 
with experiments [3, 4, 5, 6, 7, 8, 9, 10, 11, 28], quantum [16, 29], semi-classical simulations [18, 30, 31]. Surprisingly, 
even classical simulations [32, 18, 17, 22, 23, 33, 34, 35, 1] manage to reproduce such results and identify a classical 
recollision scenario as the dominant double ionization channel in the relevant intensity range. The reason for this 
£f-\ ■ success is the presence of sufficient electron-electron correlations in all such models. 

Within the dipole approximation and using soft Coulomb potentials [36, 22, 23, 33], the Hamiltonian of a one 
dimensional helium atom driven by an intense laser field is written in atomic units (a.u.) as 



H(x,y,p x ,p y ,t) = ^- + ^- + - — — - — + (x+y)E smcot, (1) 
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where x and y are respectively the position of the two electrons (along the axis of polarization) and p x and p y their 
(canonically) conjugate momenta. The laser parameters are the amplitude of the field Eq, and its frequency a>. 

Despite the two electrons are treated equally in Hamiltonian (1), a careful analysis of typical trajectories reveals 
dramatically different behaviors for each of them [35]. One electron (the outer electron) is quickly picked up by the 
field and ionized while the other one (the inner electron) is left within a competition between the laser excitation and 
the (soft) Coulomb attraction from the core. This asymmetry is explained by the dynamical organization of phase 
space by a couple of periodic orbits in the absence of the field [35]. From this competition between the laser and the 
Coulomb attraction, mainly two outcomes are possible 

• If the inner electron is sufficiently far away from the nucleus (in a region referred as the unbound area), the laser 
excitation imposes the dynamics and the electron ionizes quickly. 

• If the inner electron is sufficiently close to the nucleus (in a region referred as the bound area), it is trapped on an 
invariant torus and has no chance to ionize by itself. 



The respective size of these bound and unbound regions depend on the parameters of the laser field. To ionize an 
electron in the bound area one has to find a sufficient amount of energy to move to the unbound area. This can be 
achieved when the outer electron comes back to the core region and collides, i.e. exchanges energy, with the inner 
electron. Due to the periodicity of the laser field, the outer electron is brought back to the core nearly periodically 
since the laser is linearly polarized. Therefore the recollisions can be modeled as kicks experienced periodically by the 
inner electron. 



REDUCED MODEL FOR THE RECOLLISION DYNAMICS 

Since the dynamics in the bound area is filled by invariant tori, it is natural to perform a change of coordinates into 
action-angle variables. This (canonical) change of coordinates associates a constant action with each torus on which 
the inner electron is temporarily confined. Thus the jumps from one torus to another, experienced by the inner electron 
during the recollision process, is reduced to a change of action in the new set of coordinates. By statistical analysis 
of the recollision process together with theoretical reduced models for the inner and outer electrons, a reduced kicked 
Hamiltonian for the dynamics of the inner electron driven by recollisions was proposed [1]. Its integration yields to 
the following symplectic map : 

A„+i = A n /(l-esm(p n ), 

<Pn+i = <p„ + TV2exp(aA„ + i) + ecos<p n , 

where A n and <p„ are respectively the action and angle variables of the inner electron right before the « th kick. The 
parameter T is half of the period of the laser field, i.e. T — n/co. The parameter e, which only depends on the 
parameters of the laser field (here its intensity since we have kept the frequency fixed), simulates the variation of 
recollision efficiency as the laser intensity is varied and it is given by 

aEk 

£ = 



where the parameter a = k/ (4V2co 2 ) with k = 3.17... can be computed analytically and the parameter j3 = 

3K 2 Y 2 /(16r]L 4 co 2 ). Since there is an adjustable parameter corresponding to the interaction length between the two 
electrons (L), j3 is computed using a numerical fit (for more details, see Ref. [1]). The intensity of the laser is 
proportional to the square of the amplitude and we note that the efficiency scales proportionally to the intensity for 
low intensities, and it decreases as the inverse of its square root for larger intensities. The parameter a = — 9\/2/16 is 
obtained using a series expansion of the dynamics around the nucleus. 

Because of the finite size of the bound area, there exists a boundary invariant torus (defined as the outermost invariant 
torus analytically conjugated to a rotation (A(f),<p(f)) = (A m ,£lt + <p m ))- It defines a maximum action A m for the 
validity of the reduced model (2). Since the size of the bound area depends on the intensity, so does A m . It enables one 
to predict a ionization probability of the inner electron driven by recollisions : An electron whose dynamics is modeled 
by the map (2) is considered as ionized a soon as its action A„ becomes larger than A m after the n th kick where the 
iteration of the map (2) is stopped. 

Compared to statistical results on the double ionization probabilities obtained with the full Hamiltonian (1), 
the reduced kicked dynamics described by the map (2) agrees very well in the whole range of intensities where 
nonsequential double ionization is the dominant channel. This reduction simplifies drastically the complexity of the 
analysis : One ends up with an explicit discrete map (with two variables) instead of a continuous flow with two and a 
half degrees of freedom. In what follows we study some properties of the map (2) to connect them with the dynamics 
of the inner electron subjected to periodic collisions with the outer one. 

ANALYSIS OF THE MAP (2) 
Asymptotic behavior 

Various kicked models have been developed for atomic systems such as Rydberg atoms subjected to microwave 
fields [37, 38]. With cold atoms such kicked maps have enabled one to identify and understand accelerator modes [39]. 
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FIGURE 1. Two left panels: Finite Time Lyapunov maps of the dynamics given by Eq. (2) for laser intensities 3 x 10 W ■ cm 
(left panel) and 10 15 W- cm~ 2 (middle panel). For comparison, we also display some trajectories in the stable regions for a large 
number of iterations of the map (2). Right panel: Linear stability indicators as given by FTL exponents as a function of the laser 
field intensity. We compute the surface of stable phase space (red diamonds, left hand horizontal axis) and the mean FTL exponent 
after 10 iterations of the map (2) (blue circles, right hand horizontal axis). For comparison, we also display the variation of e with 
intensity (dotted black curve, arbitrary units). Stable initial conditions are those for which the associate FTL exponent is small at 
the end of the iteration of the map. 

For the problem at hand, because of the term T V2exp(aA„ + i) in Eq. (2), the acceleration condition is asymptotic, 
when A n — > °° (since a < 0). This asymptotic limit gives two possible values of the angle for the accelerator mode 
condition <p„ = ±7r/2 such that <p n +\ ~ <p„ and A n+ \ = A„/(l =Fe) in Eq. (2). Here £ is a positive constant and 
depending on the sign in the previous equation, A n+ \ is either larger or smaller than A„. Since the condition holds 
when A n — >• °°, only <p„ = n/2 corresponds to an asymptotic acceleration. 

This acceleration is clearly visible on a (complete) phase space of the mapping (2) and is associated with a 
burst toward very large actions around the angle n/2. However, this condition is asymptotic and thus holds only 

for \a\A„ > 1 or equivalently A„ ^> 16/ (9\/2) ~ 1 .25 which exceeds the critical action A m for ionization of the 

inner electron all over the range of considered intensities. As a consequence, accelerator modes do not contribute 
significantly for the problem at hand. They would contribute for low intensities and for rather long laser pulses. 



Linear stability 



A quick inspection of trajectories of the map (2) shows regular and chaotic areas in phase space. As the efficiency 
of the kicks (e.g., e) increases, such phase-space portraits become more sparse: elliptic (stable) islands shrink and 
the diffusion in the chaotic layer increases [1]. Using linear stability analysis such as given by finite time Lyapunov 
(FTL) exponents [40] we analyze the dynamics in the accessible phase-space for the inner electron. The way to do 
this analysis is through FTL maps which indicate the value of the exponent as a function of the initial condition in 
phase-space. Two FTL maps are depicted in Fig. 1 and compared with trajectories in the elliptic islands. FTL maps 
clearly highlight the coexistence of stable (small FTL exponent, blue regions) and chaotic (large FTL exponent, red 
regions) areas. In the chaotic layer, we also note the decrease of the FTL exponent when e decreases (by comparing 
color maps between the figures for the two laser intensities 3 x 10 14 W • cm~ 2 and 10 15 W • cirT 2 ) 

If the inner electron has initial conditions in a stable region (e.g., in one of the elliptic islands), then it will be 
impossible for the outer electron to dislodge the inner one regardless of the number of kicks (which is related to the 
pulse duration). These regions in phase space do not lead to double ionization. We investigate further the correlation 
between the (linear) stability of the map (2) and the efficiency e of the kicks, namely we compute the phase space 
surface for which the inner electron is linearly stable (small FTL exponent) and the mean value of the FTL exponent 
in the chaotic layer (large FTL exponent). All over the range of intensities where the mapping model can be used to 
simulate the collision dynamics, we see that the variation of the area of stable motion is correlated to the efficiency of 
the kicks : This area increases when the efficiency of the kicks decreases which is expected. We also note the perfect 
matching (up to a multiplicative constant) between the mean value of the FTL exponent (mpji) in the unstable region 
and the efficiency of the kicks (i.e. e <* mpji). As the FTL exponents are indicators of the (un)stability of trajectories, 
it confirms what we have already qualitatively observed looking at phase-space portraits [1]. 



CONCLUSION 

We have related linear and nonlinear properties of a two-dimensional symplectic map with the electron-electron 
collision dynamics driven by intense and short laser pulses. Accelerator modes, only accessible in the asymptotic 
limit of large actions are shown not to play a role in the dynamical process. The linear stability analysis shows the 
strong correlation between the regularity of the dynamics and the efficiency of the kicks : The stable area shrinks with 
an increase of the efficiency of the kicks, and the mean Lyapunov exponent in the unstable region is proportional to 
the efficiency of the kicks. 

ACKNOWLEDGMENTS 

We acknowledge useful discussions with P. J. Morrison and A. Kamor. CC and FM acknowledge financial support 
from the PICS program of the CNRS. This work is partially funded by NSF. 

REFERENCES 

1. F. Mauger, C. Chandre, and T. Uzer, Phys. Rev. Lett. 104, 043005 (2010); Phys. Rev. A 81, 063425 (2010). 

2. W. Becker, and H. Rottke, Contemporary Physics 49, 199-223 (2008). 

3. D. N. Fittinghoff, et al. Phys. Rev. Lett. 69, 2642-2645 (1992). 

4. K. Kondo, et al. Phys. Rev. A 48, R2531-R2533 (1993). 

5. B. Walker, et al. Phys. Rev. Lett. 73, 1227-1230 (1994). 

6. S. Larochelle, A. Talebpour, and S. L. Chin, J. Phys. B. 31, 1201-1214 (1998). 

7. T. Weber, et al. Nature 405, 658-661 (2000). 

8. C. Cornaggia, and P. Hering, Phys. Rev. A 62, 023403 (2000). 

9. C. Guo, and G. N. Gibson, Phys. Rev. A 63, 040701 (2001). 

10. M. J. DeWitt, E. Wells, and R. R. Jones, Phys. Rev. Lett. 87, 153001 (2001). 

11. J. Rudati, et al. Phys. Rev. Lett. 92, 203001 (2004). 

12. P. B. Corkum, Phys. Rev. Lett. 71, 1994-1997 (1993). 

13. K. J. Schafer, et al. Phys. Rev. Lett. 70, 1599-1602 (1993). 

14. A. Becker, and F. H. M. Faisal, J. Phys. B. 29 (1996). 

15. R. Kopold, et al. 85, 3781-3784 (2000). 

16. M. Lein, E. K. U. Gross, and V. Engel, Phys. Rev. Lett. 85, 4707-4710 (2000). 

17. K. Sacha, and B. Eckhardt, Phys. Rev. A 63, 043414 (2001). 

18. L.-B. Fu, et al. 043416 (2001). 

19. R. Panfili, and J. H. Eberly, Opt. Express 8, 431-435 (2001). 

20. I. Barna, and J. Rost, Eur. Phys. J. D 27, 287 (2003). 

21. J. Colgan, M. S. Pindzola, and F. Robicheaux, Phys. Rev. Lett. 93, 053201 (2004). 

22. P. J. Ho, et al. 94, 093002 (2005). 

23. P. J. Ho, and J. H. Eberly, Phys. Rev. Lett. 95, 193002 (2005). 

24. C. Ruiz, L. Plaja, and L. Roso, Phys. Rev. Lett. 94, 063002 (2005). 

25. D. A. Horner, et al. Phys. Rev. A 76, 030701 (2007). 

26. J. S. Prauzner-Bechcicki, et al. Phys. Rev. Lett. 98, 203002 (2007). 

27. J. Feist, et al. and J. Burgdorfer, Phys. Rev. A 11, 043420 (2008). 

28. M. Y. Ivanov, M. Spanner, and O. Smirnova, J. Mod. Opt. 52, 165-184 (2005). 

29. J. B. Watson, et al. Phys. Rev. Lett. 78, 1884-1887 (1997). 

30. J. Chen, J. H. Kim, and C. H. Nam, J. Phys. B. 36, 691 (2003). 

31. T. Brabec, M. Y. Ivanov, and P. B. Corkum, Phys. Rev. A 54, R2551-R2554 (1996). 

32. R. Panfili, and W.-C. Liu, Phys. Rev. A 67, 043402 (2003). 

33. R. Panfili, S. L. Haan, and J. H. Eberly, Phys. Rev. Lett. 89, 1 13001 (2002). 

34. J. Liu, D. F. Ye, J. Chen, and X. Liu, Phys. Rev. Lett. 99, 013003 (2007). 

35. F. Mauger, C. Chandre, and T. Uzer, Phys. Rev. Lett. 102, 173002 (2009); J. Phys. B. 42 (2009). 

36. J. Javanainen, J. H. Eberly, and Q. Su, Phys. Rev. A 38, 3430-3446 (1988). 

37. G Casati, I. Guarneri, and D. Shepelyansky, IEEE 24, 1420-1444 (1988). 

38. R. Blumel, and W. P. Reinhardt, Chaos in Atomic Physics, Cambridge U. Press, 1997. 

39. A. Buchleitner, et al. Phys. Rev. Lett. 96, 164101 (2006). 

40. P. Cvitanovic, et al. Chaos: Classical and Quantum, Niels Bohr Institute, Copenhagen, 2008, 
http : //ChaosBook . orgChaosBook .org. 



